Electromechanical instability in suspended carbon nanotubes 
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We have theoretically investigated electromechanical properties of freely suspended carbon nan- 
otubes when a current is injected into the tubes using a scanning tunneling microscope. We show 
that a shuttle-like electromechanical instability can occur if the bias voltage exceeds a dissipation- 
dependent threshold value. An instability results in large amplitude vibrations of the carbon nan- 
otube bending mode, which modify the current-voltage characteristics of the system. 



Nanoelectromechanical systems (NEMS) have been a 
subject of extensive research in recent years. The possi- 
bility of combining electrical and mechanical degrees of 
freedom on the nanoscale may give rise to technologi- 
cal advantages as well as manifestations of fundamental 
physical phenomena. From a technological point of view 
the interest is largely due to the many applications that 
may be realized using NEMS. 1 Among the many NEMS 
phenomena of considerable physical interest, we focus in 
this Letter on the "shuttle" instability that is expected 
to occur in a specific class of systems. This phenomenon 
was predicted by Gorelik et al£ and have been theoreti- 
cally investigated in a number of publications, see Shek- 
ther et al£ for a review. The typical shuttle system con- 
sists of a movable conducting grain placed in an elastic 
medium between two metallic leads. When the system is 
DC-biased over a certain threshold, an electromechanical 
instability occurs, resulting in large-amplitude mechani- 
cal vibrations and a shuttle current proportional to the 
vibration frequency. Different experimental approaches 
have been used to build and investigate various versions 
of this system^' 8 However, some problems persist and 
a true shuttle instability has not been observed. Car- 
bon nanotubes (CNT) are materials that have proved 
to be strong candidates for building NEMS. The electri- 
cal properties of CNTs depend on their detailed struc- 
ture and the mechanical properties are extraordinary. 
A very high Young's modulus and low mass give typ- 
ical resonance frequencies in the GHz-regime and the 
tubes can be mechanically deformed without breaking. 
Progress in fabrication of CNTs shows that it is possible 
to build structures, for example suspended tube a 9 ' 10 i n 
and nanorelaySfiS where CNT mechanical degrees of free- 
dom are important. The work presented here was done 
in response to recent experimental progress by LeRoy 
et al~ in measurements on suspended carbon nanotubes 
(CNT) using a scanning tunneling microscope (STM). 
Peaks in the dl/dV — V characteristics when current is 
injected into the free hanging part of the nanotube are 
attributed to inelastic tunneling of electrons that interact 
with the radial breathing mode (RBM) of the tube. This 
is a manifestation of a coupling between the electronic 
and mechanical degrees of freedom in the system, allow- 
ing for additional channels of charge transport through 



the system by means of vibronic assisted tunneling. We 
would like to point out that in addition a different elec- 
tromechanical effect — an electromechanical shuttle-like 
instability — may occur in the system. Such an instabil- 
ity occurs when the ground state is unstable with respect 
to small spatial deviations of the nanotube from its equi- 
librium position. Small displacements grow in time and 
develop into limit cycle oscillations that modify the trans- 
port properties. The shuttle instability is fundamentally 
different from vibronic assisted tunneling, in the sense 
that it results in a condensate of coherent phonons^ and 
hence large-amplitude mechanical bending mode vibra- 
tions. Below we will show that a shuttle-like instability 
may indeed occur in the suspended CNT-structure. It 
may therefore be a useful experimental system for investi- 
gating shuttling phenomena, also in the quantum regime. 
Quantum effects are expected to be important if the am- 
plitude AXq of the zero-point bending-mode oscillations 
is at least as large as the tunneling length A ~ 0.5 A - 
a length that characterizes the tunnel barrier separating 
the STM tip and the nanotube. Theoretically, the shut- 
tle system has been modeled in both the quasi-classical 4 
and the quantum mechanical limits;^ while so far the 
quantum regime has not been experimentally accessible. 
We will show below that the suspended CNT structure 
is the best candidate system for experiments that probe 
the quantum regime of shuttling. 

In the present work we investigate a system where an 
STM-tip is used to measure the current through a sus- 
pended CNT as in the experiment by LeRoy et al& We 
use the model system depicted in Figure Q and consider 
a metallic CNT (single or multi-wall) suspended over a 
trench. The nanotube is treated as a classical continuum 
beam with length /, inner radius r\ and outer radius r . 
An STM-tip is positioned near the center of the trench 
in close proximity to the nanotube; the distance from the 
undeformed tube to the tip is L ~ 1 nm. The STM acts 
as an electrode, which is connected to the tube by a tun- 
nel junction. The STM-tube separation is not fixed but 
depends on the deflection X of the tube from the equi- 
librium straight position. Both the resistance R\ and 
the capacitance C\ characterizing this junction are func- 
tions of this parameter. The ends of the tube rest on 
electrodes to which they are electrically connected via 
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Figure 1: Model system (left) and equivalent circuit (right). 
A carbon nanotube of length / is suspended over a trench; 
the tube ends are connected to metal electrodes. The tube 
can be mechanically deformed and X measures the deflection 
of the center of the tube from its equilibrium position. An 
STM-tip is positioned a distance L above the center of the 
tube and a DC-bias voltage V is applied between the STM 
and the electrodes. The CNT is treated as a metallic island 
with excess charge q and electrostatic potential ip. The STM- 
tube contact is a tunnel junction characterized by resistance 
Ri(X) and capacitance C2(X). The contacts between tube 
and electrodes are also tunnel junctions, characterized by a 
constant resistance R2 and a constant capacitance C2. 



tunnel junctions rather than ohmic contacts - the most 
likely scenario for CNT-metal junctions. The parameters 
for these junctions can be described by a single constant 
capacitance C 2 and a single constant resistance R 2 . A 
DC-bias voltage V is applied between the STM and the 
electrodes. The time scale for charge redistribution in the 
tube is fast compared to other time scales in the system 
allowing one to treat the nanotube as a metallic island 
with an excess charge q and electrostatic potential (p. A 
capacitive force F cap that depends on the potential drop 
between the STM and the CNT attracts the tube toward 
the STM-tip. 

We first discuss the possibility for an electromechan- 
ical instability to occur in our model system. For this 
purpose we restrict ourselves to the case when R\ ~ R 2 
in equilibrium. Because capacitances are not strongly 
dependent on geometry, we furthermore neglect the vari- 
ation of the capacitance with position. Since R\ depends 
exponentially on X, R2 R\ when the tube-STM sepa- 
ration decreases significantly on the scale of the tunneling 
length and R\ 3> R2 when the separation increases sig- 
nificantly. This implies that when the tube is close to the 
STM-tip, the time scale t\ for charge exchange between 
the STM and the tube is much shorter than the time 
scale T2 for charge exchange between the tube and the 
electrodes. Hence, the potential on the tube approaches 
the STM-potential. In the opposite limit, when the tube 
is far from the STM, t\ 3> t 2 the tube potential ap- 
proaches the potential of the electrodes. Consider now 
oscillations with frequency lo so that the tube position 
alternates between these limiting values. Finite values of 
the parameters lvti j2 lead to retardation, which implies 
that the potential on the tube does not have time to reach 
its equilibrium value determined solely by the tube posi- 
tion. Rather, the potential depends on the tube position 
at earlier times as well, i.e. on the tube velocity and ex- 



ceeds the equilibrium value when the tube moves away 
from the tube and vice versa. It follows that the forc e on 
the tube is correlated with the tube velocity, XF CSLp =/= 0, 
and that mechanical work will be performed on the tube 
during the oscillatory motion. 

These qualitative arguments show that the necessary 
force- velocity correlation for an instability to occur will 
develop. However, the performed work must be posi- 
tive and large enough to exceed the energy that can be 
dissipated. If this condition is met, a small deviation 
from the equilibrium position of the tube grows in time 
and results in limit cycle oscillations. The instability will 
change the current flow through the system, the nature 
of this change will be discussed below. 

To support the qualitative arguments, we have made 
quantitative calculations for our model, focussing on the 
fundamental bending mode of the doubly clamped nan- 
otube. The eigenfrequency of this mode can be obtained 
from linear elasticity theory as lo w (22A/l 2 )y/EI / pS, 
where E is the Young's modulus of the CNT, typically of 
ordeui^i^ 1 TPa , I — ? (r* — rf) is the area moment of 
inertia of the cross section, p is the mass density and S is 
the tube cross-section area. Excitation of the fundamen- 
tal bending mode can be described by a single parameter 
X(t), which denotes the deflection of the center of the 
tube from its equilibrium position. In the limit of small 
deflections (X -c I) the potential energy of the bent tube 
can be expressed as U = kX 2 /2, where k is an effec- 
tive spring constant for the tube. Hence, the equation of 
motion for the mechanical degree of freedom is 

m cB X + 1 X + kX = F cap . (1) 

The effective mass of the tube m e ff = k/uj 2 w 0.4r7i to t in 
cq ^ is proportional to the total tube mass mtot and a 
phenomenological viscous damping 7X characterized by 
a quality factor Q = Lom c tf/j has been introduced. To 
evaluate the capacitive force we use model capacitances 
for the junctions. The STM-tube junction is approxi- 
mated by a parallel plate capacitance C\ = Cq/(1—X/L), 
where Co is the typical capacitance between the STM and 
the tube and L is the equilibrium STM-tube separation. 
The tube-electrodes capacitance C2 is a constant. 

It is convenient to work with the excess charge q 
on the tube as a variable instead of the tube po- 
tential. The capacitive force then takes the form 
F cap = F (C 2 V - q) 2 /{l - aX) 2 , where F depends on 
the capacitances and a ~ 1/L. Charge is treated as 
a continuum variable, which is valid if the charge ex- 
change rate is high and the temperature not too low. 
Using an exponential form for the STM-tube resistance, 
Ri(X) = i? exp(— X/X), we derive a differential equa- 
tion for the time evolution of the charge according to 

where C S {X) = C X {X) + C 2 and C_ = d(X) - C 2 are 
the sum and difference of the junction capacitances and 
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G S (X) = l/Ri{X) + 1/R 2 and G_ = l/Ri(X) - 1/R 2 
are the sum and difference of their conductances. 

Equations ^ an d [21 define our model system, which 
we will now analyze both analytically and numerically. 
Defining (Xo,qo) as their steady state solution we lin- 
earize these equations with respect to small deviations 
from the steady state. We find that a deviation grows 
with time for fixed V if the dissipation is below a thresh- 
old value 7thr- This implies that a shuttle-like instability 
does occur in the system. However, the threshold dissi- 
pation is not a positive definite function of the system 
parameters and if it is negative no instability will occur. 
In this case the force- velocity correlation leads to intrinsic 
dissipation in the system, which suppresses any mechani- 
cal deviation and prevents an instability from developing. 
This effect originates from the position dependence of the 
STM-tube capacitance, which contributes a term to the 
force- velocity correlation that has the opposite sign com- 
pared to the term coming from the position dependence 
of the resistance. Since capacitances are not strongly de- 
pendent on geometry (L 3> A), the correlations due to 
the change in resistance usually dominates. Only in the 
limit i? 2 3> Ro does the position dependence of the ca- 
pacitance become important since the potential of the 
tube is then close to the STM-potential and the voltage 
drop over the STM-tube junction does not change appre- 
ciably with tube position. 

We can investigate the threshold dissipation analyti- 
cally in various limits. We find that the threshold dis- 
sipation goes to zero in the limits (Rq — > 0, R2 — > 0), 
(Ro — » 00, i?2 00) and Rq — > 00 while it approaches 
a negative value if R2 — > 00. First consider the limit 
(Ro — * 0,i?2 — > 0) which implies (wt\ — > 0,o;t 2 — > 0) 
and a sufficiently fast charge exchange rate for the tube 
potential to always have its equilibrium value. Re- 
tardation is negligible and no mechanical work is per- 
formed on the tube. In the limit (Rq — ► 00, R2 — ► 00), 
(u>t% — > oo,wT2 — > 00) charge exchange is a slow process 
and the tube potential is approximately constant during 
tube oscillations, making retardation a small effect also in 
this limit. Finally, in the limits Rq — ► 00 or R2 — ► 00, we 
have LUT2 3> u)t% or LOT2 3> ojti, respectively, and the po- 
tential on the tube is mainly determined by the coupling 
to the electrodes or the STM. The latter case leads to 
dissipation as discussed previously and the former leads 
to a small force-velocity correlation since corrections to 
the equilibrium potential is small. 

Our qualitative analysis suggests that optimal condi- 
tions for an instability can be found in a symmetric setup 
where Ro ~ i?2- The threshold dissipation for a typical 
set of parameters^ is plotted as a function of junction 
resistances in Figure |5J where a region of optimal con- 
ditions for the instability can be identified. This region 
corresponds to a symmetric case, where lot\ ~ LOT2 ~ 1, 
in support of the qualitative conclusion that resistances 
of the same order of magnitude are optimal. Because 
of this, we do expect poor conditions for an instability 
in the experiment by LeRoy et alS> since their junction 




Figure 2: Threshold dissipation as a function of junction re- 
sistances for V = 1 V. Dissipation is proportional to 1/Q, 
the inverse quality factor for mechanical oscillations, and an 
instability occurs if 1/Q < 1/Qthr- Large values of 1/Qthr 
correspond to optimal conditions for an instability to occur, 
and are found when the junction resistances are of the same 
order of magnitude. In the limit of highly asymmetric junc- 
tions 1/Qthr is small. 



resistances are highly asymmetric. 

The above arguments show that the occurrence of a 
shuttle-like electromechanical instability for a fixed volt- 
age crucially depends on the Q-factor of the bending 
mode vibrations and the junction parameters. Experi- 
mental Q-values of the order of 10 1 — 10 3 have been re- 
ported for CNTsjii*i£ corresponding to weak dissipation 
and an instability feasible with bias voltages smaller than 
1 V for a large range of parameters. 

In order to investigate the instability in more detail we 
solve the system of equations ^ and [3 numerically for dif- 
ferent bias voltages using a typical set of parameters. 18 
The resuts confirm the analytical conclusions, i.e. if the 
voltage exceeds a certain threshold an electromechani- 
cal instability occurs unless the threshold dissipation is 
negative. The amplitude of the limit cycle mechanical vi- 
brations is shown as a function of bias voltage in Figure 
[31 where a threshold is clearly visible when the vibration 
amplitude A becomes non-zero. The instability is associ- 
ated with different characteristics depending on the type 
of transition. The vibration amplitude is either a step 
function of voltage (hard excitation) or a smooth func- 
tion (soft excitation). In the case of a hard excitation the 
current will be a step-like function of voltage, whereas a 
soft excitation leads to a jump in the first derivative of 
the current. This can be seen in Figure 0] where the 
I — V and the dl/dV — V curves corresponding to the 
instabilities in Figure are shown. 

As we have seen, the electromechanical instability 
changes the transport properties of the system, leading 
to a step in the I — V or dl/dV — V characteristics de- 
pending on the type of transition that occurs. These 
are features that may be used to detect the instability. 
Furthermore, theoretical studies of the current noise in 
the ordinary shuttle^ have shown that the instability is 
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Figure 3: An electromechanical instability occurs when the 
bias voltage exceeds a threshold value, which results in bend- 
ing mode vibrations of the nanotube with amplitude A. The 
amplitude either increases smoothly (solid line) or in a step- 
like fashion (dashed line) as a function of voltage. Here 
Q = 50. 
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Figure 4: / — V characteristics corresponding to the soft (a) 
and the hard (b) transition in Figure [3] Insets depict the 
dl I dV — V characteristics in the vicinity of the threshold volt- 
age. The soft transition manifests itself as a jump in the value 
of dl/dV, whereas the hard transition gives rise to a peak in 
the dl/dV — V-curve. 

associated with an increase of current noised Since we 
expect a shuttle-like instability in the CNT system to re- 
sult in a qualitatively similar noise increase, noise mea- 
surements may also be used to identify the instability. A 
third approach to detecting the instability is to apply a 
small modulation to the bias voltage in analogy with the 
Sazonova et al. mixing experiment^ 

Quantum effects in shuttling have been considered 
theoreticallj*^, but they have so far not been accessible 
by experiments. With the suspended CNT system one 
can control the amplitude of the zero point oscillations 
by changing the length of the tube. It should therefore be 



possible to experimentally probe the regime where quan- 
tum effects are important by using longer tubes. An 
estimate of the ratio 2AXq/X — of order one or larger 
in the quantum regime — for different parameters shows 
it to be in the range (0.5 — 1) (///im) 1 / 2 for a single wall 
carbon nanotube. This is a larger value than can be esti- 
mated for other realizations of the shuttle e.g. one finds 
2AX /X ~ 0.2 - 0.3 in the C 60 -SET device used by Park 
et al& The suspended CNT-shuttle is so far the system 
where quantum corrections will be most significant. 

The typical STM-tube separation is of order 1 nm and, 
consequently, surface forces between the tube and the 
STM will be important. However, the influence of these 
forces is mainly quantitative, it does not change the qual- 
itative conclusion that an instability occurs. The scale 
of these forces for interactions between a CNT and a 
metallic electrode istf of the order of a few nN. In the 
suspended CNT-structure the scale of the elastic forces 
can be changed by orders of magnitude by varying the 
tube length and thickness. Consequently, the importance 
of the surface forces relative to the elastic forces depends 
on the tube geometry. 

To conclude, we propose that a shuttle-like electrome- 
chanical instability can be experimentally studied by us- 
ing suspended carbon nanotubes in tunneling contact 
with supporting metal electrodes and an STM tip. We 
support this claim with an analysis of a classical model 
of such a system and seek the optimal conditions for an 
instability to occur. The shuttle-like instability is most 
likely to occur in the symmetric case, when all the tun- 
neling resistances are of the same order of magnitude. 
We demonstrate that the effect of the instability is to 
change the I — V characteristics and, finally, we show 
that the suspended nanotube system is the first, where 
the quantum regime may be approached experimentally. 

Acknowledgment. Financial support from the 
Swedish Research Council and the Swedish Foundation 
for Strategic Research is gratefully acknowledged. This 
work was in part supported by EC FP6 funding (con- 
tract no. FP6-2004-IST-003673, CANEL). This publica- 
tion reflects the views of the authors and not necessarily 
those of the EC. The Community is not liable for any use 
that may be made of the information contained herein. 



* Electronic address: mjonsson@fy.chalmers.se 

1 Cleland, A. N. Foundations of Nanomechanics, Springer- 
Verlag: Berlin, 2003. 

2 Gorelik, L. Y.; Isacsson, A.; Voinova, M. V. Kasemo, B.; 
Shekther, R. I.; Jonson, M. Phys. Rev. Lett. 1998, 80, 
4526. 

3 Shekther, R. I.; Gorelik, L. Y.; Galperin, Yu.; Isacsson, A.; 
Jonson, M.; J. Phys. Cond. Mat., 2003, 15, R441. 

4 Fedorets, D.; Gorelik, L. Y.; Shekther, R. I.; Jonson, M. 
Europhys. Lett. 2002, 58, 99. 

5 Fedorets, D.; Gorelik, L. Y.; Shekther, R. I.; Jonson, M. 



Phys. Rev. Lett. 2004, 92, 166801. 

6 Novotny, T.; Donarini, A.; Jauho, A. P. Phys. Rev. Lett. 

2003, 90, 256801. 

7 Scheible, D. V.; Blick, R. Appl. Phys. Lett. 2004, 84, 4632. 

8 Park, H.; Park, J.; Lim, A. K. L.; Andersson, E. H.; 
Alivisatos, A. P.; McEuen, P. M. Nature 2000, 407, 57. 

9 LeRoy, B. J.; Lemay, S. G.; Kong, J.; Dekker, C. Nature 

2004, 432, 371. 

Lee, S. W.; Lee, D. S.; Yu, H. Y.; Campbell, E. E. B.; 
Park, Y. W. Appl. Phys. A 2004, 78, 283. 

1 V. Sazonova, Y. Yaish, H. Ustimel, D. Roundy, T. A. Arias, 



■5 



and P. M. McEuen, Nature 431, 284 (2004). 
Lee, S. W.; Lee, D. S.; Morjan, R. E.; Jhang, S. H. Sven- 
ingsson, M; Nerushev, O. A.; Park, Y. W.; Campbell, E. 
E. B.; Nano Lett. 2004, 4, 2027. 

Poncharal, P.; Wang, Z. L.; Ugarte, D.; de Heer, W. A. 
Science 1999, 283, 1513. 

Wong, E. W.; Sheehan, P. E.; Lieber, C. M. Science 1997 
277, 1971. 



Treacy, M. M. J.; Ebbesen, T. W.; Gibson, J. M. Nature 
1996, 381, 678. 

Isacsson, A.; Nord, T. Europhys. Lett. 2004, 66, 708. 
Jonsson, L. M.; Nord, T.; Viefers, S.; Kinaret, J. M. J. 
Appl. Phys. 2004, 96, 629. 

Co = C 2 = 5 • 10~ 18 F, / = 200 nm, r D = 5 nm, and n = 1 
nm. 



